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ABSTRACT 

As part of a long-term project to revisit the kinematics and dynamics of the large 
disc galaxies of the Local Group, we present the first deep, wide-held (~ 42 7 x 56 7 ) 
3D-spectroscopic survey of the ionized gas disc of Messier 33. Fabry-Perot interfer¬ 
ometry has been used to map its Ha distribution and kinematics at unprecedented 
angular resolution (<3 77 ) and resolving power (~ 12600), with the 1.6 m telescope 
at the Observatoire du Mont Megantic. The ionized gas distribution follows a com¬ 
plex, large-scale spiral structure, unsurprisingly coincident with the already-known 
spiral structures of the neutral and molecular gas discs. The kinematical analysis of 
the velocity held shows that the rotation center of the Ha disc is distant from the 
photometric center by ~ 170 pc (sky projected distance) and that the kinematical 
major-axis position angle and disc inclination are in excellent agreement with photo¬ 
metric values. The Ha rotation curve agrees very well with the Hi rotation curves for 
0 < R < 6.5 kpc, but the Ha velocities are 10 — 20 km s _1 higher for R > 6.5 kpc. 
The reason for this discrepancy is not well understood. The velocity dispersion prohle 
is relatively hat around 16 km s _1 , which is at the low end of velocity dispersions 
of nearby star-forming galactic discs. A strong relation is also found between the Ha 
velocity dispersion and the Ha intensity. Mass models were obtained using the Ha 
rotation curve but, as expected, the dark matter halo’s parameters are not very well 
constrained since the optical rotation curve only extends out to 8 kpc. 

Key words: techniques: interferometric, integral held, Fabry-Perot - galaxies: indi¬ 
vidual: M33 - galaxies: kinematics and dynamics - galaxies: ISM: Bubbles, HII regions 
- cosmology: dark matter 


INTRODUCTION 


Late-type spiral and dwarf irregular (dlrr) galaxies present 
an extended Hi emission. The Hi emission usually extends 
to larger radii than the optical emission. This makes of Hi 
a powerful tool for probing the dark matter (DM) domi¬ 
nated regions and for characterizing the parameters of the 
flat part of the rotation curves (RCs). However, it is the ris¬ 
ing inner part of the RCs that constrains the parameters of 
the mass models (see e.g. Blais-Ouellette et al. 19991. Two- 
dimensional Ha observations are ideally suited to give more 
accurate velocities in the inner parts of galaxies having a 


better spatial resolution than the Hi (Amram et al. 1992 
|1994| |Swaters et al.||2000| ). The ideal approach is then to 
combine high spatial resolution Ha observations in the inner 
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parts and sensitive lower spatial resolution Hi observations 
in the outer parts to derive the most accurate mass mod¬ 
els for the luminous disk and the dark halo. The objective 
of this project is to perform a new mass distribution model 
of Messier 33 (M33), combining high-sensitivity Ha and HI 
interferometric data. In particular, this article presents the 
first Ha survey devoted to the large-scale distribution and 
kinematics of the M33 ionized gas disk. 


M33 is the second most luminous spiral (SA(s)cd) 
galaxy in our neighborhood after M31. With an absolute 


magnitude of My~-19.3, it presents several arms (Sandage 


& Humphreys 19801. Although two main arms are well 


known, M33 has not a clearly defined grand-design pattern. 
A system of seven spiral arms is pointed out by |Ivanov H] 


Kunchev (1985). In the optical bands, it could be consid¬ 


ered as a flocculent spiral galaxy, but UV and IR observa¬ 
tions show more prominent arms. The M33 profile presents 
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Table 1. Optical parameters of M 33. 


Parameter 

Value 

Source 

Morphological type 

SA(s)cd 

RC3 

R.A. (2000) 

01 h 33 m 33.1 s 

RC3 

Dec. (2000) 

+30° 39' 18" 

RC3 

Optical inclination, i 

52°±3° 

ww 

Optical major axis, PA 

22.5°±1° 

ww 

Apparent magnitude, my 

5.28 

RC3 

Absolute magnitude, My 

-19.34 


(J-K) 

0.89 

JTH 

Optical radius, R 25 ( ; ) 

35.4 ± 1.0 

RC3 

Systemic velocity, V sys 

-179 ±3 km s - 1 

RC3 


References:WW: |Warner et al. | (|1973[l , RC3:|de Vaucouleurs 
et al.| (|1991[) , JTH: |Jarrett et al.| (|2003[> 


Table 2. Distances for M 33. 


Method 

(indicator) 

Distance Modulus 
m-M [mag] 

Distance 

[Mpc] 

Source 

TRGB 

24.54 ± 0.06 

0.81 

(i) 

TRGB 

24.75 ± 0.10 

0.89 

( 2 ) 

TRGB 

24.64 ± 0.15 

0.85 

(3) 

PNLF 

24.62 ± 0.25 

0.84 

(4) 

Cepheids 

24.64 ± 0.06 

0.85 

(5) 

Cepheids 

24.70 ± 0.13 

0.87 

( 6 ) 

Cepheids 

24.52 ± 0.14 

0.80 

(7) 

Cepheids 

24.55 ± 0.28 

0.81 

( 8 ) 

Cepheids 

24.58 ± 0.10 

0.84 

(9) 

Cepheids 

24.62 ± 0.01 

0.84 

( 10 ) 

Mean adopted distance 

0.84 ±0.04 
1" = 4 pc 



References: (1) |McConnachie et al.|( 

2005), (2) |Stonkute et al.| 

(2008|), (3) Galleti et al. (2004), (4) 

Magrini et al.||2000|, (5) 

Saha et al. (2006 

1 , ( 6 ) Paturel et al. 

(|2002), (7) |McConnachie | 

et al. (2004J, ( 8 ) 

An et al. (|2007|), (9) (Freedman et al.||2001|, 


(10)|Gieren et ah] (|2013[l 


a small bulge-like component in the IR bands. Its bulge and 


nucleus are the subject of numerous studies (e.g.: 

McLean & 

Liu 1996;|Lauer et al. 1998; Gordon et al. 1999; 

Corbelli & 

Walterbos 

2007 

1. The optical parameters of M33 are given 

in Table [1 




In the literature, one can find distances to M33 varying 
from 700 to 940 kpc. Table [2] gives recent distance mea¬ 
surements for M33. The adopted distance for this study 
is D = 0.84 Mpc (scale 1’ = 244 pc). It is the mean dis¬ 
tance estimated by using resolved sources techniques such as 
Cepheids, Planetary Nebulae Luminosity Function (PNLF) 
and Tip of the Red Giant Branch (TRGB). The optical disk 
of M33 has a scale length of ~9.2 / (2.25 kpc) in the V band 
(Kent 1987| [Guidoni et al.||198l| ) and its optical radius is 
R 25 = 35.4 (RC3). A cutoff in the radial surface brightness 
profile appears at ~ SO 7 (~ 9 kpc) in the I band (Ferguson 
|et al.|2007| >. 

Mapping the environment of the Local Group galaxies, 
as is done by the “Pan-Andromeda Archaeological Survey” 
(PAndAS) ( |McConnachie et al.||2009|) , confirms that there 
was many mergers and interactions between them. The dis¬ 
covery of dwarf galaxies around the Milky Way and M31 and 
the tidal streams between M31 and M33 (PAnDAS) confirm 
our ideas about galaxy formation. The particular structure 
(star streams) seen in M33 could be associated with this 
history of mergers and interactions. In fact, many of the 


structures presented in McConnachie et al. (2009 2010); 


Cockcroft et al. (20131 and Wolfe et al. (20131 around M31 


and M33 can only be explained by these interactions. Deep 
Ha observations of M33 reveal the presence of low density 
Hu regions outside the optical disk (Hoopes et al. 2001). 


This suggests recent star formation activity, possibly due to 
recent interactions. 

Studying the kinematics of such a galaxy will provide 
a better understanding of the contribution of dark matter 
and of the best functional form describing its distribution. 
Still today, the cusp-core problem remains as one compares 
observations to predictions, especially for dwarf systems and 
this, despite the numerous studies on the DM distribution in 
galaxies in the last 30 years. While, on large scales, N-body 
cosmological simulations reproduce well the observations, it 
is more problematic at galaxy scales. The NFW profile, de¬ 
rived from those simulations, predicts a cuspy distribution in 


the central parts of the DM halos (e.g. 

Fukushige & Makino 

1997; Moore et al. 1999; Navarro et al. 

1997, 2010, Ishiyama 

et al.||2013 

l, while observations especially of dwarf systems 


show more a core distribution (Oh et al. 2011). Those re¬ 


sults can be explained by the gravitational potential related 
to the gas in those simulations, since the gas, which is impor¬ 
tant in the inner parts, is not accurately introduced in those 
simulations. Phenomena such as stellar feedbacks, galactic 
winds or massive clumps are not often well handled and re¬ 


produced ([Goerdt et al.|2010| Inoue & Saitoh 2011| Ogiya & 

Mori 2011 Pontzen & Governato|2012 Teyssier et al.|2013). 

Several studies of the kinematics of M33 exist (e.g. 

Cor- 

belli & Salucci|2000 Corbelli 2003, Putman et al. 2009J 

Cor- 

belli et al.|2014 

i, mainly based on Hi. It appears from those 


studies that M33 is dark matter dominated (the dark mat¬ 
ter mass is > 10 times larger than the mass of the stellar 
disk) and that its Hi disk is strongly warped. The gas in the 


disk of M33 is estimated ~ 1.4-3x10 M@ (Corbelli 


2003 


|Putm an et al.|2009| ). In the mass model derived by |Corbelli| 


120031, the stellar mass is estimated at ~ 3 —6x 10' Mq and 
the dark halo mass at ~5xl0 10 Mq. Most of the M33 kine- 
matical studies use Hi or CO (H 2 ), with a probable beam 
smearing impact on the innermost parts of the RC due to the 
lower spatial resolution of the Hi data. Among the few Ha 
(optic al) studies, the results presented by |Carranza et al.| 


(1968) show small velocity dispersions from 5 km s -x to 
9 km s' 1 in the disk. However, they argue that these values 
could even be lower if taking into account the instrumental 
corrections. The mass-to-light ratio is an important param¬ 
eter in the determination of the dark matter halo shape. In 


M33, Boulesteix & Monnet (1970) show very low values of 
M/L in the central (< 50') part of the galaxy and increasing 
in the outer parts of the disk. The kinematics of the inner 


regions was studied using the warm gas (Ha) (Boulesteix & 
Monnet 1970 Boulesteix et al. 1974| Zaritsky et al. 1989I 


and the cold molecular gas (H 2 via CO) (Wilson & Scoville 
1989[ IGratier et al.|2010[ |2012| |Kramer et al.||2013| |Druard 


et al. 2014). 


The precise determination of the rising parts of RCs 
with high resolution data and a better estimate of the M/L 
ratio for the stellar disk could define more accurately the 
shape of its dark matter halo. High resolution optical ob¬ 
servations of a nearby galaxy such as M33 is complex in 
view of the large size of the galaxy compared to the small 
field of view of high resolution integral field spectroscopic 
instruments. For M33, there is a lack of high resolution op- 
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tical data available for the kinematical study. An exception 
is the study by Corbelli & Walterbos (2007) on the central 
5' using the gas and stellar kinematics obtained by long-slit 
spectroscopy. 

Mass models are sensitive to the rising part of the RCs 


(see, e.g Amram et al. 1996 Swaters et al. 1999 Blais 
Ouellette et al.||1999l. High resolution Ha RCs give a more 


accurate determination of the kinematical parameters for 
the derivation of the RCs and subsequently a more accurate 
determination of the mass model parameters . The high res¬ 
olution Fabry-Perot (FP) integral field observations at Ha 
with a resolution > 10000 give optimal kinematical data 

for the optical disk. The Ha RC can be used to test the 
shape of the DM halo allowing us to compare the derived 
DM distribution with CDM predictions. 

To address these problems, this study presents Fabry- 
Perot (FP) mapping of M33 obtained at the Observatoire du 
mont Megantic (OMM). R elafio et al.| p013) have studied 
the Spectral Energy Distribution (SED) of the Hu regions 
of M33 and the star formation rate (SFR) and star formation 


efficiency (SFE) have been investigated by Cardan et al. 

(20071 and |Kramer et ah (20111. More than 1000 Hll regions 
can be resolved by the Ha observations; |Courtes et ah| ( |1987] ) 
gave a catalogue of 748 Hll regions. Observing those regions 
allows us to derive the ionised gas (optical) kinematics of 
M33. Determination of the M33 kinematics with a spatial 
resolution <3” (~ 12 pc) using the Ha velocity field and 
the derivation of an accurate RC in the inner parts are the 
main goals of this paper. The 3D data will be used to derive 
mass models with a dark halo component (ISO and NFW). 

Section [2] describes the data acquisition and the instru¬ 
mentation used while section |3] discusses the data reduction 
techniques. Section [4] details the kinematical analysis and 
section[5]gives the details of the mass models analysis. A dis¬ 
cussion of the M33 velocity dispersion, a comparison with 
other studies and the results from the mass modeling can 
be found in Section [6] Finally, a summary and the general 
conclusions are given in section [7]. 


2 FABRY-PEROT OBSERVATIONS 

2.1 Telescope and instrument configuration 

The observations took place at the 1.6-m telescope of 
the Observatoire du mont Megantic (OMM, Quebec), in 
September 2012. A scanning Fabry-Perot (FP) etalon inter¬ 
ferometer has been used during the observations with the 
device IXON888, a commercial Andor EM-CCD camera of 
1024x1024 pixels . The details of this camera, based on e2v 
chips, are given in Table [3] 

The IXON888 camera provides a large field-of-view 
(FOV) and was set to Electron Multiplying (EM) mode, 
14 bits read-out resolution and its detector was cooled to 
—85 K during the observations. The camera clocks, gains 
and read-out speeds were adjusted in order to reduce the 
noise. 

The order-sorter filter is a narrow band interference fil¬ 
ter, centered at A c = 6557A (nearly at the systemic velocity 
~ —180 km s -1 ) with a FWHM of 30A. Its maximum trans¬ 
mission is ~80%. The interference order of the FP interfer¬ 
ometer is p = 765 at Ha. The FP has a Free Spectral Range 


Table 3. Photon counting cameras at OMM. 


Cameras 

IXON888 “ 

Pixels size 

0.84"x 0.84" 

Active Pixels 

1024x1024 

Quantum Efficiency (QE) 

> 90% 

Cooling 

-85 °C 

Max. Readout Speed 

10 MHz 

RON 

< 1 e~ with EM gain 

Detectors 

CCD201-20 6 

CIC C levels 

5 - 8 x 10~ 3d 


1 = 4 pc, 0.8 = 3.3 pc at a distance of 0.84 Mpc 

“ http://www.andor.com 

b http://www.e2v.com/products-and-services 
Clock- induc e d cha rges per pixel per frame 


Daigle et al. 


1 2009 I 


RON : Read-Out Noise given per pixel per frame 


(FSR) of 8.16A (~373 km s _1 ), which has been scanned 
through 48 channels, corresponding to a spectral sampling 
of ~ 0.17A (7.8 km s -1 ). The finesse / of the Fabry-Perot 
etalon provides the spectral resolution AA = FSR/f. Our 
observations were done with a mean finesse of / = 16.5, as 
determined from datacubes of a Neon calibration lamp. The 
FWHM spectral resolution is thus AX max ~0.53A (resolv¬ 
ing power of ~ 12620 at Ha). This corresponds to a FWHM 
instrumental broadening of 23 km s -1 (dispersion of ~ 10 
km s^ 1 ) at the scanning wavelength of the observation of 
6558.8 A. 


2.2 Data acquisition 

The wide field-of-view of ~14 / x 14 ; allows to map the bright 
inner disk regions of Messier 33 with only ten pointings. The 
different pointings overlap by a few arcseconds to allow the 


derivation of one final large mosaic with no gaps (S 3.31. Fig¬ 


ure [I] shows those pointings, whose centre coordinates are 
listed in Table [4] Notice that the central field has been ob¬ 
served more than once to yield very deep Ha data for the in¬ 
nermost regions of Messier 33. A region of the sky, free from 
apparent Hll regions from Messier 33, has been observed 
before and after each galaxy pointing whose duration was 
larger than 30 minutes. Those acquisitions were done with 
the same interference filter as the one used for the galaxy, in 
order to perform the subtraction of the night-sky emission 
lines. 

The data were obtained by operating the camera at 
two seconds exposure per frame. A gap of 0.4s between two 
consecutive channels was necessary to move the reflective 
plates of the FP etalon in order to avoid overlapping frames 
during the fast transfer. 

Dark current, gain, flat field and noise calibration ob¬ 
servations have been acquired at the beginning and the end 
of each night. The “Dark” observation consisted in a series of 
at least fifty images during which the detector is exposed for 
two seconds without light, and the gain is set to its largest 
amplitude, like during the observations of every sky/galaxy 
fields-of-view. The “Gain” observations consisted of about 
200 frames, again acquired with the largest gain value. The 
“Noise” observations consisted in an integration of zero sec¬ 
ond, with the lowest gain value. A minimum of fifty frames 
were collected in order to calculate the read out noise of the 
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Right Ascension (J2000) 


Figure 1. Mosaic of the 10 Fabry-Perot pointings in the direction of M33, shown with the sky field to the NW. The FOV of each field 
is ~ 1 4' x 14\ with an overlap of ~ 5". The background image is a R band archive image from KPNO i]Hoopes et al.|2001). 


CCDs. The noise, gain and dark were observed in off-light 
mode in the dome. 

The total exposure time performed for each field is listed 
in Table [4] In total, the time spent to integrate on the 
Messier 33 fields was >30 hours. The seeing of the observa¬ 
tions was < 3”. To perform the wavelength calibration, the 
system (filter+lenses+FP+camera) has been illuminated by 
a Neon (Ne) lamp. A narrow band filter centered on 6598A 
and of FWHM=16.3A was used to isolate the Ne line at 
A=6598.95A. 


3 DATA REDUCTION 

3.1 Wavelength calibration, spectral smoothing 

The raw frames of the observations contain interferograms 
giving the information on the number of photons per frame 
per channel and per cycle. We used Interactive Data Lan¬ 
guage (IDL) routines to integrate raw 2D files into a 3D dat- 
acube. A phase correction consisting in shifting every pixels 
such that they are at the same wavelength across the field 
has been applied to the raw datacubes to yield wavelength- 
calibrated datacubes. For that purpose, a phase map has 
been derived from a Ne line calibration observation. Dat¬ 
acubes are then wavelength-sorted and corrected for guid¬ 
ing shifts and cosmic rays. The ’’Noise”, ’’Dark” and ’’Gain” 


observations (12.2 1 were then used to correct and calibrate 
the detector response. 

A Hanning filter with a width of three channels has then 
been applied to the spectra to increase the sensitivity. This 
reduction step is the same as the one applied to the Virgo 


Cluster, GHASP or SINGS samples ( 

Chemin et al. 

2006 

Epinat et al. 2008 ?| Daigle et al.!;2006 

Dicaire et al. 

to 

o 

O' 

00 


The same process is used for the sky observations, whose 
datacubes have been subtracted of all the other fields to 
produce M33-only Ha? emission line datacubes. For pointings 
having more than one observation, all night-sky corrected 
and wavelength calibrated datacubes have been combined 
into a single datacube, to increase the sensitvity. 


3.2 Airmass and ghosts corrections 


Because the center of the interference rings of the obser¬ 
vations is near the center of the detector, there are some 
reflections about the optical axis, called ghosts. The size 
and the intensity of the ghost depend on the shape of the 
region which produces the effect. We have used the method 
described in Epinat (2009) to reduce the number of pixels 


dominated by those ghosts. The typical residuals after ghost 
removal are about 10% of the initial ghosts flux. 

Then, an airmass correction has to be applied. Indeed, 
differences in airmass during the multiple sessions of obser- 
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Table 4. Characteristics of the observations 


Fields ID 

RA 

DEC 

T 

S 

am 

v h 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

M33_01a 

1:33:59 

+30:40:48.54 

2.0 

2.7 

1.1 

19.6 

M33_01b 

1:33:59 

+30:40:48.54 

1.4 

2.2 

1.3 

17.5 

M33_01c 

1:33:59 

+30:40:48.54 

2.0 

2.2 

1.1 

17.5 

M33_02a 

1:33:22 

+30:29:22.27 

1.4 

2.8 

1.1 

17.0 

M33_02b 

1:33:22 

+30:29:22.27 

2.7 

2.5 

1.3 

17.0 

M33-03 

1:34:29 

+30:51:52.41 

2.4 

2.4 

1.3 

19.3 

M33_04 

1:33:06 

+30:40:16.96 

2.5 

2.4 

1.1 

19.2 

M33-05 

1:34:50 

+30:40:30.72 

2.5 

2.3 

1.4 

18.9 

M33_06 

1:34:15 

+30:28:52.98 

3.0 

2.3 

1.1 

19.0 

M33-07 

1:33:35 

+30:52:04.25 

2.6 

2.1 

1.3 

18.5 

M33_08 

1:32:33 

+30:28:16.19 

2.7 

2.1 

1.1 

18.5 

M33-09 

1:33:26 

+30:18:26.23 

2.4 

2.2 

1.1 

17.4 

M33 10 

1:34:18 

+31:02:01.11 

2.6 

2.5 

1.1 

17.2 

Sky 

1:32:14 

+31:00:09.48 

0.2 



17.5 

Dark 



0.1 




Noise 



0.1 




Gain 



0.1 




Flat 

Dome 

Dome 

0.1 





Columns notes: (1) Identification of the field. (2) and (3) 
Coordinates of the field center (J2000). Each field-of-view is 
14' x 14'. Observations presented in this paper have been taken 
in September 2012. (4) Exposure time in hours, (5) Seeing 
FWHM during observations. (6) Air mass correction factor. (7) 
Vh is the heliocentric velocity correction (km s —1 ). 


vations can affect the number of counts received by the de¬ 
tector, which implies field-to-field sensitivity variations. In 
Table |4j we give the airmass correction factor (AM) which 
has to be applied to each field in order to get the counts 
equivalent to an airmass of 1.0. Starting from the observed 
counts Cobs at a given airmass, the true value of the counts 
Ctrue is given by: 

Ctrue = 10“°- 4AX X C obs (1) 


The parameter A represents the Ha extinction coefficient. 
Our Fabry-Perot system has already been used by|Hlavacek-| 


Larrondo (20091 to determine the value of that coefficient. 


We used their value of 1.03. The X quantity is given by 
X = AM —1.0, with the airmass correction factors AM listed 
in Tab. □ 


3.3 Mosaicing and binning the datacubes 

Due to the large angular size of Messier 33, there are two 
possibilities to generate the final set of moment maps (Ha 
integrated emission, line-of-sight velocity and velocity dis¬ 
persion maps, continuum map). A first one is to compute for 
each of the ten pointings its own moment maps, and com¬ 
bine all the maps to produce a single field-of-view moment 
map. The second approach is first to combine all datacubes 
into a mosaic datacube, and then generate moment maps. 
The first approach is easier to do, but it introduces more 
errors in radial velocities in regions of overlap. Instead, we 
think it is better to combine and correct the Ha spectra di¬ 
rectly, than radial velocities or velocity dispersions. We thus 
chose the second approach, which is a more powerful pro¬ 
cess, though less straightforward to implement. Combining 
the data cubes increases the signal to noise per pixel, which 


results in increasing the numbers of pixels for which we can 
measure a line. To minimize the errors in the overlap regions, 
astrometry information have been added in each single cube 
headers. The different steps in the process are : 

(i) White light image (sum of all channels) 

(ii) Accurate astrometry on the white light image 

(iii) Field positioning from astrometry 

(iv) Generation of an exposure map with the sum of ex¬ 
posures 

(v) Accurate spectral calibration, taking into account he¬ 
liocentric correction 

(vi) Generation for each channel of a mosaic 

(vii) Weighting of each channel by the exposure map 
(viii) Combination of all channels to generate the mosaic 

cube 


In order to process projections and make an accurate 
mosaic, some parts of the Montage packages from IPAC (In¬ 
frared Processing and Analysis Center) have been used to 
develop our tools. 

Once the mosaic datacube has been obtained, an adap¬ 
tive spatial binning has been performed. That adaptive spa¬ 


tial binning is based on Voronoi 3D binning (e.g. Daigle et al. 
2006 Chemin et al.||2006| Dicaire et al.|2008 l. The Voronoi 


technique consists in combining pixels to larger bins up to 
a given value of S/N (signal-to-noise ratio) or more. The 
pixels/bins with a S/N lower than the threshold S/N are 
combined with the neighbouring pixels until the threshold 
S/N value is reached. A S/N of 7 has been targeted for this 
study. This particularly allows to preserve the angular reso¬ 
lution where the S/N is high. 

Figure [2] gives the channel maps of the data cube used 
to determine the kinematic al parameters of M33. The detec¬ 
tions plotted in Figure [2] are those above 3 sigmas. It rep¬ 
resents the wavelength variation from 6556.5A to 6560.6A 
(from channel 11 to channel 35 in our data) where the chan¬ 
nels width is 0.18 A. The Hu regions appear progressively 
from the North East side to the South-West side. 


3.4 Moment maps derivation 

The integrated emission field is the 0-th moment map de¬ 
rived from the spectra and the velocity field the first mo¬ 
ment map. Radial velocities are given in the heliocentric rest 
frame. Making the assumption that the PSF has a gaussian 
profile, the velocity dispersion field is the second moment 
map: 


& cor r 


E i A, 2 F(A i ) 

EiF(Ai) 


A 2 °PSF 


( 2 ) 


where F(Ai) is the flux corrected for the continuum level 
at wavelength A i, A is the barycenter of the emission line, 
and (7 psf is the instrumental dispersion. Velocity disper¬ 
sions have not been corrected for thermal broadening of the 
medium, nor for the natural width of the Ha line. While the 
natural width of the Ha line remains negligible (~3 km s _ ), 
the reason for not correcting for the thermal dispersion of 
the gas was that the temperature of the ionized interstel¬ 
lar medium of Messier 33 is not known accurately. A typical 
broadening of ~ 10 km s _1 is often given in the literature, for 
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Figure 2. Ha channel maps. The radial velocity (km s ) is indicated in the upper left corner of each channel and the cross indicate 
the kinematical center. 


a temperature of ~ 10 4 K for the ionized gas. However, we 
have observed many bins of dispersion (after instrumental 
broadening correction) lower than this usual value of ~ 10 
km s _1 , implying that a typical temperature of ~ 10 4 K 
cannot be valid everywhere. 

Finally, the moment maps have been cleaned to get rid 
of all possible unrealistic patterns, like those from reflection, 
noise and background emission residuals. For that purpose, 
we have masked all pixels with data values lower than twenty 


counts. We are left with maps having 546941 independent 
bins. 

3.5 Flux calibration 

The flux calibration that turns Fabry-Perot data counts into 
surface brightness (in erg/s/cm 2 /arcsec 2 ) is given by: 

SB = Cst x SBfp 


( 3 ) 
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Right Ascension (J2000) 


Figure 3. Selected Hu regions for the Ha flux calibration over¬ 
plotted on the H a line flux map obtained from our set of obser¬ 
vations. The aperture size around each region is the same as in 
[Relano et ah (|2013|>■ 


Here, Cst is a calibration constant, SB is the corresponding 
surface brightness value, and SBf p is the calibrated value 
of FP count /pixel/s. The FP fluxes calibration using a Ha? 
map is a linear relation. We have used the catalog of Hu 
regions of Relano et ah](j2013) as reference fluxes to calibrate 
our interferometric data. Figure[3]shows the selected regions 
and their associated aperture size used to integrate the Ha 
counts. For each region the total integrated flux is computed 
using the counts and the mean exposure time. 

The calibration constant for the M33 Ha FP im¬ 
age is given by 1 count/pixel/s = 2.45 ± 0.03 x 

lCU 17 erg/s/cm 2 /arcsec 2 . The typical error that can be no¬ 
ticed on the flux determination at an aperture used, is 
~2.4 x 10 _lb erg/s/cm 2 /arcsec 2 . It takes into account the 
errors on the Ha flux provided in Relano et al. (2013) and 
the dispersion of the fit. Figure [4] shows the comparison bew- 
teen our Fabry-Perot calibrated data, and the |Re!ano et al.] 
12013) data in units of flux, erg/s/cm 2 . 


4 Ha DISTRIBUTION AND KINEMATICS 
4.1 Ha maps 

The observations obtained at the OMM give 13 cubes on 
ten fields. Those cubes were used to build a mosaic cube 



Calibrated flux (1CT 12 x erg/s/cm 2 ) 


Figure 4. Comparison between the derived Fabry-Perot Ha 
fluxes and the reference Ha fluxes of |Relano et al.| { j201 used 
for the FP flux calibration. Square symbols represent data from 
individual Hu regions (Fig.|3j and the solid line is the linear fit 
to the data. 
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Figure 5. Ha map of M33 derived from our Fabry-Perot data. 
The CO peak contours are in red (120 K km/s) from Tosaki] 
|et al.| ([2011) and the HI peak contours are in blue (9500 and 
1200 K km/s) from [Chemin et al.| ( [2012| ) . The black sketchs are 
drawn to follow the Hii regions and trace the multi arms profile 
of M33. These arms are close to those shown in |Li &; Henning| 
(|201 1[> or |Boulesteix et al~] (|1974[). 


as described in section |3.3| Figure [5] shows the Ha emis¬ 
sion of M33 in the ~42 ; x 56 / field. The contours are from 
Hi and CO emission. We can see that while the Ha emis¬ 
sion follows roughly the Hi emission, it is even more so for 
the CO emission. The Hi and CO structures seem to fol¬ 
low the arms described by the Hu regions. The final cube 
was produced with the same FSR and same spectral reso¬ 
lution as the small cubes observed. The Ha monochromatic 
image (moment zero map) is presented at the top right of 
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Figure [6] and can be compared to the WISE I NIR image 
on the left. The discrete Hll regions have different sizes and 
shapes, filled and clear shell regions as described by |Relano| 


et al. (20131. Two main strong arms are clearly defined along 


with the multi arms structure as presented by [Boulesteix] 
etTaE] ( fl974| >. 

The velocity field (first moment map) at the bottom left 
of Figure [6] was obtained by using the zeroth moment map 
as a mask. The bins of the voronoi binning with data values 
greater than 20 counts are shown. This criteria was chosen 
in order to avoid all probable unreal pattern introduced by 
ghosts of strong Ha regions, noise and/or background. The 
velocity dispersion map (second moment) is shown in the 
bottom right of Figure [6] 

For this study, the data will be used mainly to derive 
the overall kinematics and derive the optical rotation curve. 
Figure [?] shows that the resolution of this data set (~ 3”) 
could be used to study the detailed kinematics of Hu re¬ 
gions, shells, cavities, bubbles, filaments, loops, outflows and 
ring-like structures. This will be done in another publication. 
Except within such local features, the large-scale Ha veloc¬ 
ity field of M33 seems regular, typical of nearby disk galax¬ 
ies, without significant apparent twist of the major axis, and 


with mild streamings inherent to spiral arms (see sect. 4.2.41. 

Figure [^provides a zoom of the ~ 8 ; x 8' central regions 
of M33. When compared to the WISE I NIR image to the 
lop left, it can be seen clearly that emission is detected all 
the way to the very center of the galaxy. The velocity dis¬ 
persion map shows that brighter HII regions exhibit larger 
dispersion (20-30 km s^ 1 ), though the largest dispersions 
(up to 60 km/s) are only seen in bins with the faintest Ha 
emission. The map also shows that the inner regions exhibit 
smaller dispersion (see Section 6 for more details). 


4.2 Kinematics 

We performed a kinematical analysis in order to recover the 
set of kinematical parameters and the best possible rotation 
curve with meaningful uncertainties. 


4-2.1 Initial parameters 

The ROTCUR algorithm implemented in the reduction 


package GIPSY (Groningen Image Processing System; Vo- 
|gelaar fc~T erlouw 200T| has first been used to find the mean 
initial values of the kinematical parameters, and to ver¬ 
ify that no significant warp of the ionized gas disc exists. 
ROTCUR is based on the tilted-ring model described in 
Begeman (1987). Starting from initial values of i = 51° and 
PA = 200° , V 3 ys and the rotation center (xo, yo) are first de¬ 
termined. A second run allowed us to derive the inclination 
and PA profiles. From those profiles, we measured an incli¬ 
nation i = 52° ± 2°, and a PA = 202° ±2°, both values being 
very close to the optical morphological values (see Table [l]). 
In both cases, the errors are determined using the devia¬ 
tion around the mean values. No trend has been detected 
in the inclination and position angle profiles, and the small 
standard deviations are clear indications that no significant 
warp of the Ha disk exists (at least in the regions covered 
by our observations), as is usually the case for the optical 
disks (warps are mainly seen in the extended Hi disks). 


4-2.2 Rotating disk model 

Since no warp is detected in the Ha disc using ROTCUR 
(tilted ring model), we used a model in which the gas is 
supposed to lie in a unique plane in order to derive the kine¬ 
matical parameters and the rotation curve. Thus we use the 
whole 2D information to derive the projection parameters 
and their uncertainties. The model is explained in |Epinat| 
et al. (2008) and was used for the GHASP sample data 


analysis. When the vertical motion velocities are not con¬ 
sidered, the observed line-of-sight velocities are expressed 
as: 

Vobs = V sys -I- (Vrot(R) cos(0) + Vexp(R) sin(0)) sin(i); (4) 

where R and 9 are the polar coordinates in the plane of the 
galaxy, i is the disk inclination, V ro t(R) is the azimuthal 
velocity (i.e. the rotation curve) and V e x P (R) is the radial 
velocity in the galaxy plane (often referred to as the expan¬ 
sion velocity). Defining the kinematical position angle PA as 
the anticlockwise angle between the North and the direction 
of the receding side, the azimuthal angle 9 can be deduced 


at each position (more details in the annexes of Epinat et al. 

|2008| ). 

With the hypothesis that the radial velocities of the 
ionized gas are negligible with respect to the azimuthal 
rotation, the observed velocities become V 0 b s = V ays + 
Vrot(R) cos(0) sin(i). We therefore build a 2D model with 
a set of projection parameters (center, position angle and 
inclination) and a set of kinematics parameters describing 
the rotation curve and the motion of the galaxy with re¬ 
spect to the Earth (systemic velocity). All these parameters 
have no dependency with radius, contrary to what is done 
in tilted ring models (e.g. Begeman 1987). 

The rotation curve we used in our model is described 


by the Zhao function (Kravtsov et al. 19981 with reduced 
parameters: 


V rot = V t 


(R/rt) E 
1 + (R/rt) a 


(5) 


where r t and Vt define the transition (“turnover”) radius and 
velocity, a and g describe the sharpness of this transition and 
the shape before and after the transition. Therefore, the 2D 
model is described by a set of 9 parameters (i, PA, X cen , 
Keen, V sys , Vt, rt, g, a). The optimization starts with the 
initial parameters i, PA, Vsys and the rotation center pre¬ 
viously found using ROTCUR, which are compatible with 
morphological parameters. The method uses a \ minimiza¬ 
tion calling the IDL LMFIT routine based on the Levenberg- 
Marquardt method in order to find the best fit model. The 
uncertainties on the parameters are derived using a Monte 
Carlo method based on the power spectrum of the residual 
velocity map (see Epinat et al. 2008 for details). 


4-2.3 Derived kinematical parameters 

With thousands of degrees of freedom, the optimized model 
converges rapidly towards a stable solution. The optimized 
center is found at R..A. = 01 h 33 m 54.1 B , DEC. = 2>Q 0 2)QA2". 
That location is ~ 42” (168 pc) from the photometric cen¬ 
tre (sky projected distance), towards the NE direction, with 
an angle of 60° with respect to the semi-major axis of the 
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Figure 7. Same as in Figure[6] zoomed on the star forming region located at a=01h32m43s, <5=30°24’57”. 


approaching side. This offset corresponds to a deprojected 
distance of 63” (252 pc) in the galaxy plane of Messier 33. 
As a comparison, sky-projected offsets between photometric 
and kinematical centres of bright spirals in the Virgo cluster 


of galaxies have been found between 70 and 800 pc (Chemin 
|et a l.] [2006j ). The offset we find for Messier 33 is thus com¬ 
parable with other spirals, at the low end of the distribution 
for other galaxies. 


The derived systemic velocity is V Bys = —178 ± 3 
km s _1 , the position angle PA = 202 ±4° and the disc incli¬ 
nation i = 52 ± 2°. These parameters are in excellent agree¬ 
ment with those found with ROTCUR, and with literature 
values. The values for the other parameters are g = 1.2±0.7 
and a — 1.01 ± 0.6; those parameters being in good agree¬ 
ment with best fits parameters range of rotating discs de¬ 
scribed by Kravtsov et al. (1998). Rotation can be approxi¬ 
mated to a linear function of radius in the inner parts where 
R < rt = l.^and velocities V < Vt = 50 km s^ 1 . 


4-2-4 Rotation curve 

Rotation curves for the entire disk, and for the separate ap¬ 
proaching and receding sides of the disk have been extracted 
using the projection parameters obtained from the model fit¬ 
ting described above. As detailed in Epinat et al. (2008]), the 
rotation curve is extracted in rings whose width is optimized 
by gathering a minimum of fifteen velocity measurements 
in the velocity field within each ring. Because the velocity 
field contains more than 500000 independent pixels, it was 
very easy to define rings with such a number of pixels. As 
a result, the rotation curves were incredibly resolved, with 
thousands of rings for 0 < R < 35 ; . Therefore, for practi¬ 
cal reasons we decided to rebin the resulting curves at radii 
regularly spaced by 5 , and of 5 width. The top panel of 
Figure [lO] shows the Ha rotation curves for the approaching 
and receding sides, while the bottom panel shows the global 
rotation curve, as fitted using both sides of the disk simulta¬ 
neously. The adopted rotation curve is given in Table[5]with 
the associated velocity errors AV. The adopted rotation ve¬ 
locity uncertainty is given by AV = y/e 2 + |V a — V r | 2 /4. 
The term e is the dispersion around the mean value within 
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Figure 8. M33 central ~8 / field - same caption as in Figure [5] The datapoints presented are the barycenter after an adaptive binning 
with S/N=7. The cross indicate the kinematical center. 


Table 5. Sample of the Ha rotation curve of M33. Full table available online. 


Rad 

(i) 

V 

(2) 

AV 

(3) 

(7 

(4) 

Act 

(5) 

Rad 

(6) 

V 

(7) 

AV 

(8) 

(7 

0) 

Act 

(10) 

Rad 

(ii) 

V 

(12) 

AV 

(13) 

(7 

(14) 

Act 

(15) 

0.08 

6 

4 

14 

10 

11.08 

87 

4 

16 

5 

22.08 

104 

5 

17 

7 

0.17 

9 

4 

12 

6 

11.17 

86 

3 

16 

5 

22.17 

105 

6 

17 

7 

0.25 

13 

3 

12 

4 

11.25 

86 

3 

16 

6 

22.25 

106 

8 

18 

7 

0.33 

18 

5 

14 

5 

11.33 

86 

3 

16 

6 

22.33 

107 

7 

19 

8 

0.42 

21 

4 

15 

4 

11.42 

86 

4 

15 

6 

22.42 

106 

7 

18 

8 

0.50 

23 

4 

16 

4 

11.50 

85 

3 

15 

6 

22.50 

105 

5 

17 

7 

0.58 

25 

3 

15 

4 

11.58 

83 

3 

15 

6 

22.58 

107 

5 

16 

7 

0.67 

26 

3 

16 

5 

11.67 

85 

3 

16 

6 

22.67 

107 

6 

17 

8 

0.75 

31 

3 

16 

5 

11.75 

86 

3 

16 

6 

22.75 

105 

7 

16 

7 

0.83 

32 

3 

15 

5 

11.83 

85 

3 

16 

6 

22.83 

107 

5 

15 

7 

0.92 

35 

3 

15 

4 

11.92 

87 

3 

16 

6 

22.92 

108 

5 

16 

7 

1.00 

35 

3 

15 

4 

12.00 

89 

3 

17 

7 

23.00 

109 

5 

15 

6 

2.00 

49 

4 

14 

4 

13.00 

94 

3 

18 

8 

24.00 

110 

5 

17 

8 

3.00 

56 

3 

14 

5 

14.00 

96 

4 

16 

6 

25.00 

114 

4 

16 

7 

4.00 

61 

3 

15 

5 

15.00 

98 

7 

15 

7 

26.00 

108 

19 

17 

7 

5.00 

65 

6 

16 

5 

16.00 

97 

5 

16 

6 

27.00 

110 

14 

14 

8 

6.00 

72 

5 

17 

5 

17.00 

99 

4 

17 

6 

28.00 

130 

3 

19 

12 

7.00 

78 

4 

17 

6 

18.00 

99 

3 

17 

7 

29.00 

125 

5 

19 

9 

8.00 

80 

3 

17 

6 

19.00 

107 

3 

17 

7 

30.00 

121 

3 

14 

10 

9.00 

81 

4 

17 

6 

20.00 

104 

5 

16 

6 

31.00 

122 

4 

15 

11 

10.0 

81 

4 

15 

6 

21.00 

108 

3 

15 

8 

32.00 

128 

4 

22 

21 


Notes : Column (1): radius in arcmin, (2) the rotation velocities in km s , (3) errors on V (4): u ( velocity dispersion) profile and (5) 
errors on the velocity dispersion. The following columns (up to 15) have the same definitions as the previous columns 
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Figure 9. Same as in Fig. zoomed on the star forming region NGC 604 


each ring (statistical error for the both sides velocity cal¬ 
culation). The second part |V a — V r | /2 is the systematic 
uncertainty that expresses the asymmetry between rotation 
velocities for the approaching (V a ) and receding (V r ) disk 
halves. The formal statistical error is usually smaller than 
the systematic error. The RC used for the mass models is the 
adopted RC of Figure |To| (bottom panel). Only the points 
where data are present on both sides are used. As shown 
in Figure [To] (top panel), the RC on the approaching side 
goes out to ^ 8.5 kpc and only out to ~ 8.1 kpc on the ap¬ 
proaching side. So, the adopted RC stops at ~ 8.1 kpc. The 
RC was smoothed at a binning of 5 to get constant step in 
the mass modeling. The Ha rotation curve exhibits a regu¬ 
lar rising gradient within the inner R = 8 kpc, reaching a 
maximum velocity of ~ 125 km s -1 at the last data points. 
Many wiggles are obviously seen, as probable consequences 
of the crossing of the spiral arms of Messier 33 or due to 
co-rotation effects. The axisymmetry of the rotation is very 
good. Indeed, the most significant differences (~ 35 km s -1 ) 
between V a and V r are only observed in a narrow range 
around R = 6.5 kpc. The axisymmetric velocity field model, 
as deduced from the adopted rotation curve with all kine- 


matical parameters kept constant with radius, and its sub¬ 
traction to the observed velocity field are shown in Fig. E3 
The distribution of residuals is centered on 0 km s , with 
a standard deviation of 8 km s _ , implying a very accurate 
kinematical model for most of the disk. We note that locally, 
bright Ha emission in spiral arms can exhibit larger velocity 
differences. For instance, the star forming region at R.A. = 
01 h 33 m 15.4 s , DEC. = 30°56 , 48” has an average residual 
of ~ 10 km s" . This likely shows the limit of the axisym¬ 
metric rotation model at those angular resolution (5”, i.e 
20 pc). Our observations are indeed sensitive to very local, 
non-axisymmetric motions inherent to such star forming re¬ 
gions (expanding shell, etc.), which cannot be modeled cor¬ 
rectly by the large-scale rotation of the galaxy. Notice also 
that larger residuals in spiral arms indicate the presence of 
asymmetric or streaming motions. Such motions are often 
observed in disk galaxies, but cannot be modeled by the ax¬ 
isymmetric rotation. Finally, larger residuals at the outskirts 
of the Ha disk could indicate that some Hu regions may not 
necessarily lie in the main equatorial plane of inclination i 
= 52°, so that deprojection of their velocities has not been 
performed correctly. 
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Figure 10. Ha rotation curve of Messier 33. The top panel shows the curves for the approaching and receding disk sides (blue downward 
triangles and red upward triangles, respectively). The adopted rotation curve of M33 is shown in the bottom panel. 



Rodius (orcmin) 


Figure 12. Ha velocity dispersion profile of M33. Dotted lines 
display values at + /- one standard deviation derived in each radial 
bin. The bottom panel shows the number of independent bins 
used to calculate the dispersion. 


4-2.5 Velocity dispersions 

The velocity dispersion profile has been derived by az- 
imuthally averaging the veloci ty d ispersions er in annuli of 
5 width, in 5 steps (Figure 12 I. Annuli have been com¬ 
puted using the projection parameters determined from the 


kinematical model (sect. 4.2.3 1 The a profile is almost con¬ 
sistent with a flat profile. The mean dispersion is 16 km s^ 1 . 
Like for the rotation curve, wiggles are detected. They are 
likely caused by increased dispersion when crossing the spi¬ 
ral arms. They are nonetheless not significant to make a 
noticeable increase of the scatter. We measure a standard 
deviation of the disperion profile of 2 km s _ . The pro¬ 
file marginally decreases from R=1.5 kpc towards the cen¬ 
tre. From R ~ 7.5 kpc, the velocity dispersions increase to 
~ 20—25 km s -1 . This radius is the location of the beginning 
of the warp of the Hi disk where the twist of the position 
angle starts ( Corbelli fe Salucci|2000 1. However, it would be 
interesting to have Ha data with good SNR beyond these 
radius to better understand the Ha dispersions behaviour in 
regions where the Hi disk warp is more pronounced. 


5 MASS MODELING 

The Ha rotation curve describes accurately the velocity gra¬ 
dient in the center of galaxies, usually better than any other 
kinematical tracer. Such optical high-resolution is crucial to 
test different inner shapes of dark matter haloes, like cuspy 
or shallow models. This section only focuses on the mod¬ 
elling of the mass distribution of Messier 33 within the inner 
8 kpc from our newly derived Ha rotation curve. However, 
with a RC only derived out to 8 kpc, we do not expect 
strong constraints on the halo’s parameters. We postpone 
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Figure 11. The left image gives the axisymmetric model obtained with the kinematical parameters obtained with the Zhao model (See 
Sect. |Tg. 2 ( . The map was masked with the Ha monochromatic map (top right map in figure[6|. The right image gives the residual map 
(observed-model). 


to a forthcoming article a more complete modelling of the 
mass distribution from a more extended rotation curve that 
will merge our inner Ha rotation curve with a new Hi rota¬ 
tion curve for the outer regions of Messier 33 (Kam et al., 
in preparation). 


5.1 Luminous mass components 

5.1.1 The neutral gaseous disc 


Provisional data on the Hi gas component have been pre¬ 
sented in [Chemin et al. (20121. Those observations will be 
fully described in a future Kam et al.’s paper. Figure [13] 
presents (left-hand panel) the Hi mass density profile over¬ 
laid on the ionized gas brightness profile (in arbitrary units). 
Both profiles have been derived with the task ELLINT in 
GIPSY. The Hi disc mass is ~ 2 x 10 9 Mpj. The surface 
density only slightly varies around 8 Mg pc 
kpc, then decreases at larger radii. 


0 - 

-2 


within R = 7 


5.1.2 The bulge-disk decomposition 

For the stellar contribution, the surface brightness profile is 
derived from the Spitzer/IRAC 3.6 pm data. The archive mo¬ 
saic file from guaranteed time observations of Robert Gehrz, 
Observer Program ID 5 was used (Gehrz & Willner||20041. 


After removing the stars and the background, the CHAN¬ 
DRA’S CIACrj tools have been used to derive the profile of 
Figure [l3| The Spitzer/IRAC 3.6 pm profile is compared to 
the WISE I (3.4 pm) data. The profile presented in the right- 
hand panel of Figure [l3| (red), reaches a surface brightness of 
P 3.6 ~ 21.0 mag.arcsec -2 at a radius of 21 arcmin (~5 kpc), 
which is similar to the WISE I profile (blue). The surface 
brightness p cor , corrected for the inclination and dust, is 
obtained by: 

Mcor = Mobs + A a + 2.5 x b A log(cos(i)) (6) 


where the correction parameters A\ and b A are tabulated in 


Graham & Worley (2008). The mean inclination of 52° is 


used for the density correction. In this correction, the influ¬ 
ence of PAH at 3.3 pm is considered as very weak and the 
maximum dust correction used in the IR comes from the J 
band ( |Graham fe Worley 20081. In the right panel of Fig¬ 
ure m the central 3 arcmin shows a small spheroidal compo¬ 
nent in the 3.6 pm band. This small central bulge has been 
discussed by many authors (Regan & Vogel 1994; Gebhardt 


et al. 2001 Minniti et al. 1994 Seigar 20111. Kent (19871 


found that the nucleus in M33 is similar to a point source 
and the rising inner 3* part of the surface brightness profile 


1 http://cxc.harvard.edu/ciao/ 
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Figure 13. Left panel: Hi surface density profile (blue line) and Ha brightness profile (red line, in arbitrary units). Right panel: in red 
the M33 surface brightness profile using the Spitzer/IRAC 3.6/im and in blue using the WISE I (Tom Jarrett, private communication) 
data. The dashed-dotted line is the bulge component, the dashed line is the disk component, the other line is the contribution of both 
components. 


suggests the presence of a small bulge. [Cor bclli & S alucci| 
12007) found for M33 a bulge component extending up to 
1.7 kpc and Seigarj |201 1| only to 0.39 kpc. In view of those 
different results, it is clear that a disk-bulge decomposition 
has to be done using the IR profile from this study. 

The best fit of the bulge-disk decomposition, shown in 
Figure |l3l is obtained using an exponential disk with a Sersic 
model for the bulge. The black dot-dot-dot-dash line gives 
the best fit of the sum of the decomposition. The disk is 
described by: 


With these values, our study leads to a bulge-to-disk ratio 
of B/D=0.04 ± 0.012 that is in agreement with B/D=0.03 
obtained by |Seigar| (2011) with IRAC data at 3.6 /im. The 
bulge component is subtracted from the total profile in order 
to get the disk contribution. 

Stars in a disk have a vertical thickness. For the ver¬ 
tical distribution of t he ste llar com ponent, w e adopted a 
sech 2 (z/za ) law ( van der Kruit fe Scarle 1081 ■. We used a 
vertical scale height of ~ 365 pc, which is ~20% of the stellar 
disk scale length. 


Mr) = IM> + 1-10857 (7) 

tid 

where ho is the central surface brightness and Rd the scale- 
length of the disk. The disk parameters found are Rd = 
1.82 ±0.02 kpc (slightly smaller than the scale-length in the 
optical) and ho = 18.01 ± 0.03 mag.arcsec -2 . As seen in 
Figure |13[ the Wise I profile has a slightly shorter scale- 
length. The disk parameters found for the Wise I profile are 
Rd = 1.70 kpc and ho = 17.60 mag.arcsec -2 . 

The bulge is described by : 


/x(r) 


fie + 2.5 b n 



( 8 ) 


where, /x e is the effective surface brightness at R e , the effec¬ 
tive radius. R e defines the radius that contains half of the 
total light. The parameter n determines the “curvature” of 
the luminosity profile. b n is defined as b n . = 1.9992n — 0.3271 
for 0.5 < n < 10 (CapaccioTi 1 989| ). The best fit in Figure [l3| 
gives He = 20.8 ± 0.3 mag.arcsec -2 , n = 1.12 ± 0.1 and an 
effective radius R e = 0.35 ± 0.05 kpc for the IRAC profile 
and fie = 21.5 ± 0.4 mag.arcsec -2 , n = 1.20 ± 0.05 and an 
effective radius R e = 0.40 ± 0.05 kpc for the Wise profile. 


5.1.3 The mass-to-light ratio 


With such a small bulge-to-disk ratio (and consequently Re), 
M33 can be considered as nearly a pure disk galaxy. In this 
paper, both cases are considered; the pure disk case and the 
disk ± small bulge case. However, with such a small bulge, 
they should be fairly similar. The color mass-to-light ratios 
(M/L) are defined separately for the bulge and the disk and 
are used to obtain the actual total stellar mass contribution 


(9) 


using the method described in Oh et al. (2008): 

T 3.6 =B 3.6 x T iC + ^3.6 

where Y*' 6 is the M/L at 3.6/xm, Yfy the M/L in the K 
band and B 3 ' 6 and A 3 ' 6 the correction coefficients. For the 
K band M/L, we used the relation log(T^) = 1.46(J — K) — 
1.38 taken from de Blok et al. ( 2008|), ob tained by using an 
extrapolation of Bell & de Jong (20011, where the stellar 


mass synthesis uses a Salpeter initial mass functio n (IMF). 
The T*' 6 for nearby galaxies is obtained using Y(/ (de Blok 
et al. [2008): 

T 3 - 6 = 0.92 x Tf - 0.05 (10) 

The K band M/L is obtained using the 2MASS (J-K) 
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color computed by Jarrett et al. (2003). The (J-K) color has 
been computed separately for the bulge in the inner part 
and for the disk in the outer parts. The first data points 
(J-K) are used for the bulge. For the pure disk case, the 
mean color from 65 to 570 arcsec has been used. The (J-K) 
decomposition gives for the bulge and disk 0.94 ±0.03 and 
0.86 ±0.03 respectively. However, it is quite likely that the 
bulge color is underestimated, being contaminated by disk 
light. Those values give the mean M/L in the 3.6 p,m band : 
Y))' 6 = 0.72 ±0.1 for the disk and Y 3 ' 6 = 0.80 ±0.1 for the 


bulge. The effective mass density profile is obtained by (Oh 
et al. 2008): 


E[M 0 pc“ 2 ] = Y 3 ' 6 x 10 


-0.4X(/1 3 .6-24.8). 


( 11 ) 


where YJ' 6 is the mass-to-light ratio in the Spitzer/IRAC 
3.6 micron band. The density profile is used in GIPSY to 
compute the contribution of the stellar component. 


5.2 Dark matter halo density profile 


The total rotation velocity is given by: 


V r 2 ot =V 2 ±V 2 ±v£m 


( 12 ) 


where V? is the contribution of the stars, V 2 „ the contribu- 
* 2 E 
tion of the gas component and V DM the contribution of the 

dark matter halo. The dark matter contribution is required 
to explain the outermost flat part of the rotation curves in 
galaxies (Bosma 1978 Carignan & Freeman 19851. The dark 
matter distribution can be dehned by different types of den¬ 
sity profile. We will limit this study to the most commonly 
used halo density profiles, the pseudo-isothermal (ISO) and 
the Navarro, Frenk and White (NFW) halo distributions, 
which show the largest differences in the inner parts, where 
the RC is well defined by the Ha data. 


5.2.1 ISO density profile 

The pseudo-isothermal (ISO) dark matter halo is a core- 
dominated type of halo. The ISO density profile is given by: 

piso{R) = 1 + P °_Rj 2 ( 13 ) 

where po is the central density and R c the core radius of the 
halo. The velocity contribution of a ISO halo is given by: 

Viso(R) = y / 47rGp 0 R?(l-^atan(|-)) (14) 


5.2.2 NFW density profile 

The NFW model is derived from A CDM simulations 


(Navarro et al. 

1996 1997) 

"universal halo 

) is known e 


R 1 law (de B lok|2010 ) in the innermost regions. The NFW 
halo density profile is described by: 


Pnfw{R ) = _ g _ P * _^. 2 (15) 

r s r s 1 


Table 6. Ha mass model results with DM halos for the pure disk 
case. 


Halo Model Params 

Best Fit 

Color M/L* 

(i) 

(2) 

(3) 

(4) 

ISO 

P0 

6.75 ±0.57 

7.46 ± 0.43 


Rc 

24.23 ± 17.02 

11.13 ±0.40 


Y 

0.81 ± 0.01 

0.72 


Xr 

0.88 

1.16 

NFW 

R200 

240.00 ± 188.34 

300.26 ± 106.55 


c 

4.41 ± 4.15 

2.37 ± 0.91 


Y 

0.48 ±0.11 

0.72 


Xr 

0.87 

1.56 


*M/L fixed by the color of the disk; 

po, the central DM density, is given in units of 10~ 3 M 0 /pc 3 ; 
R c and R 200 are in kpc. 


where pi ~ 3 Hq / ( 8nG ) is the critical density for closure of 
the universe and Rs is a scale radius. The velocity contri¬ 
bution corresponding to this halo is given by: 


Vnfw(R-) — V200 


ln(l ± cx) — cx/(1 ± cx) 
x(ln(l ± c) - c/(l ± c))] 


(16) 


where V200 is the velocity at the virial radius R200, c = 
R200/Rs gives the concentration parameter of the halo and 
x is defined as R/Rs. The relation betwen V200 and R200 is 
given by: 


V200 

R200 x Ho 

100 

(17) 

the 

Hubble constant 

taken as 


where Ho i 
Ho = 72 km s _1 Mpc _1 ( Hinshaw et al.|2009 l. 


5.2.3 Results 


The Ha rotation curve of Figure 10 (bottom panel) will be 


used for the mass modeling. Figure [14] shows the models 
using the ISO (top) and the NFW (bottom) DM distribu¬ 
tions for the pure disk case and Figure [l5| for the bulge-disk 
decomposition. The left panels of the figures give the best 
fit models and the right panels, the models with the mass- 
to-light ratio constrained by the IR. color and population 
synthesis models pT et al. 2008). Both fits use Levenberg- 
Marquardt least-squares fitting techniques. At the bottom 
of each model, the mean residuals (observation - model) are 
represented by a black line with the same error bars as the 
velocity errors of the top panels. The part colored in pink 
gives the dispersion of the residuals around the black regres¬ 
sion line. 

Tableland Table[7]give the results of the mass models: 
column (1) gives the type of halo density profile used; col¬ 
umn (2) gives the parameters of the halo. In the tables, Y 
gives the mass-to-light ratios obtained from the fits and \ r 
the goodness of the fit. Column (3) shows the results using 
the best fits and column(4) the results when the mass-to- 
light ratios are kept fixed at the value obtained using the 
(J-K) color and population synthesis models. In Table [ 7 ] 
the parameter for the bulge mass-to-light ratio is given by 
T 5 . The results of the mass models will be discussed in the 
next section. 
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Figure 14. M33 mass models with a DM halo for the pure disk case. Models on the left are best fits while models on the right have 
a fixed M/L determined by the IR color (see sec. ED- The top panels represent the models with the ISO and the bottom panels those 
with the NFW density distributions. The stellar disk contribution is in broken red, the Hi contribution in blue, the halo contribution 
in broken green and the sum of the 3 components is shown as a continuous red line. At the bottom of each model the mean residuals 
(observation - model) are presented. 


6 DISCUSSION 


6.1 M33 velocity dispersion 


Our new dataset is ideal to study the internal velocity dis¬ 
persion of HII regions in M33, and the mean velocity dis- 
perion of the galaxy. In Figure |16| we compare the mean 
velocity dispersion of M33 to a subsample of 151 nearby 


star-forming galacitc disks from the GHASP survey ( Epinat 
et al.|2008 l studied in Epinat et al. (20101 It is obvious how 
the ionized interstellar medium of Messier 33 appears dy¬ 
namically colder (lower velocity dispersion) than in other 
galaxies. We have divided the sample into three classes of 
size: small galaxies with R 25 <7.5 kpc, intermediate size 
galaxies with 7.5 < R 25 < 11.5 kpc and large galaxies with 
R 25 > 11.5 kpc. Though being within the class of interme¬ 
diate galaxy size, Messier 33 curiously behaves like smaller 
galaxies that have lower mass and absolute magnitude. 

Giant extragalactic Hu regions seen in gas-rich spiral 


and dlrr galaxies are regions of strong star formation with 
a size ranging between 0.1 and 1 kpc, a Ho luminosity of 
10 39 — 10 4l) erg/s, a mean density of 1 — 10 cm -3 and an 
ionized mass of 10 4 — 10 J M@, which embeds a population 
of 100 — 200 ionizing stars (|Kennicutt et al.Jl984, Kcrmicutt 


1984). In addition, the velocity dispersion of the gas presents 


distinct kinematics with subsonic and supersonic emission¬ 
line widths ([Smith & Weedman 1970,4971). Uniform expan¬ 
sion of Hu regions into a medium of constant density results 
in subsonic expansion velocities (Spitzer & Tomasko 19681. 
Supersonic motions are inferred for velocity dispersion typi¬ 
cally larger than (kT/m) 0 ' 5 ~ 10 — 13 km s _1 , where k is the 
Boltzmann constant, m the typical hydrogen mass weighted 
by the molecular content and T = 10 4 K, the characteristic 
temperature of the Hu region. Supersonic velocity disper¬ 
sions ranging from 15 to 40 km s _ are observed and give 


rise to different interpretations (see discussion in Bordalo 


& Telles (20101. In the "champagne" model (Tenorio-Tagle 
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Figure 15. Same as in Figure [l4| but with the bulge contribution (in cyan) 
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Figure 16. Mean velocity dispersion of Messier 33 (star) and of other nearby star-forming disc galaxies from [Epinat et al.| ( |2010| . On 
the right, the mass is given by M = V^ iax R 25 /G, where V m ax is the maximum velocity of each Ho; rotation curve. 



































Kinematics and Mass Modeling of MSS: Ha Observations 


19 


Table 7. Ho; mass model results with DM halos, using the bulge- 
disk decomposition. 


Halo Model Params 

Best Fit 

Color M/L* 

(i) 

(2) 

(3) 

(4) 

ISO 

PO 

6.41 ± 0.60 

12.52 ±0.53 


Rc 

119.50± > 200 

6.32 ±0.40 


r d 

0.82 ±0.01 

0.72 


r b 

1.32 ±0.06 

0.80 


Xr 

0.85 

1.34 

NFW 

R200 

211.55 ±0.72 

2.78 ±0.78 


c 

8.67 ± 0.12 

0.76 ±0.45 


T d 

0.43 ±0.01 

0.72 


Tf, 

1.01 ± 0.11 

0.80 


Xr 

0.75 

1.43 



350 

280 

210 w 
$ 

5 

140 z 

70 

0 


*M/L fixed by the color of the disk and of the bulge; 

po, the central DM density is given in units of 10 —3 Mg/pc 3 ; 

R c and R 200 are in kpc. 


Figure 17. Corrected velocity dispersion (<r 3 or = — ap SF ), 

computed using a single Gaussian, vs integrated line intensity for 
the whole field-of-view of M33. The number of pixels having a 
given intensity and dispersion is given by the color scale. 


1979), the expanding Hu bubble bursts thought the sur¬ 


face of the surrounding molecular cloud where the sharp gas 
density discontinuity generates a shock wave that acceler¬ 
ates the gas to supersonic speeds. Density gradient within 
the HII regions may generate supersonic motions prior to 
their acceleration through the surface of the cloud (|Mazurek| 
1982 1. 

Our data confirm that the region NGC 604 (Figure [ 9 ]), 
like many other giant Hu regions, displays both single Gaus¬ 
sian and complex line profiles ( Munoz-Tunon et al.|[l996| 
Yang et al.|1996 l. Single Gaussian profiles essentially come 
from the bright regions and show the broadening mechanics 
due to self-gravitation. More complex lines come from faint 
regions and emanate from the wind-driven mechanical en¬ 
ergy injection due to massive stars that produce expanding 
shells, cavities and bubbles, filaments and outflows, loops 
and ring-like regions. 

A velocity dispersion versus intensity (a — I) diagram 
provides diagnostics used to separate the main broadening 
mechanisms detected in the emission lines of giant Hu re¬ 
gions, such as differentiating the broadening produced by 
virial motions resulting from the total gravitational poten¬ 
tial and that resulting from the superposition of shells and 
loops generated by massive stars, which ends up dispersing 
the parent clouds (|Munoz-Tunon et al. 1996 |Yang et al. 


1996 |Martfnez-Delgado et al.||2007| |Moiseev fc Lozinskaya 


20121. Such a diagram is given for the whole disk of the 


galaxy in Figure [17] We used the total intensity in the line 
rather than the intensity peak because the total intensity 
is independent from the spectral resolution (which is not 
the case of the intensity peak). In the following, we use in¬ 
tensity to refer to the total line intensity. Figure [lT] shows 
the strong relation between the Ha velocity dispersion and 
the Ha intensity. Only low intensity regions show a broad 
range of velocity dispersions. In other words, the range of 
velocity dispersions decreases when intensity increases. Re¬ 
gions showing larger velocity dispersions arise from diffuse 
emitting areas rather than from intense Hn regions. Sub¬ 
sonic dispersions are only observed in low intensity regions, 
whereas supersonic motions are both observed in the low 
surface brightness medium and the brightest star forming 


regions. These latter have a roughly constant velocity disper¬ 
sion (20-30 km s _1 ), while the largest supersonic dispersions 
(40-60 km s" ) are only seen among the lowest intensities. 


As described in|Munoz-Tunon et al. 

|1996|), in the frame 

of the CSM for stellar cluster formation 

Tenorio-Tagle et al. 


19931, the gravitational collapse fragments the gas clouds 


and first forms the low-mass stars. The bow shocks and 
wakes caused by their stellar winds suspend the collapse of 
the cloud and communicate the stellar velocity dispersion to 
the surrounding gas. The core of the cloud is thus virialized 
and the supersonic gas velocity dispersion traces its gravita¬ 
tional potential, despite dissipation. Massive stars form later 
in a cloud at equilibrium in which the velocity dispersion of 
the gas is constant and does not depend on the intensity of 
the newly formed giant Hu regions. In the a — I diagram, 
this corresponds to the horizontal area (with the almost con¬ 
stant velocity dispersion of 20±3 km s _ ) displaying a broad 
range of intensities. 

As massive stars evolve, strong mechanical energy 
sweeps the ISM into shells displaying supersonic velocities 
higher than the ambient stellar velocity dispersion. The ve¬ 
locity dispersion of these shells decreases with age while their 
luminosity increases: young and faint shells produce a large 
range of velocity dispersions and low intensities while older 
and bright shells display a narrower velocity dispersion am¬ 
plitude and higher intensities. If the shells are embedded in 
the virialized core, their velocities cannot be lower than the 
stellar velocity dispersion but, if the shells are outside the 
core, their velocities can decrease to lower values, down to 
subsonic velocities. In the a — I diagram, these shell phases, 
ages and locations are found within the broad range of veloc¬ 
ity dispersions at low intensities. To fully understand how 
this last area of the a — I diagram is filled by expanding 
shells, this zone should be understood as the result of dif¬ 
ferent projection effects. 

On the one hand, emission from a line-of-sight passing 
across the center of an expanding spherical shell of gas will 
present the largest velocity difference across the shell due 
to the large distance separating its two edges and an inter¬ 
mediate intensity. On the other hand, a line-of-sight going 
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Figure 18. Comparison of Ha and Hi rotation curves. Blue rep¬ 
resents the Hi (squares) RC from |Corbelli et al.| ( |2014| ) ; black is 
for|Gratier et aiT7| (|201Q[ ) VLA data; green circles are for [ Corbelli 
|fe Schneider| ( |1997[ ) data; green crosses for |Newton] ( | 1980[ ) data 
and red symbols are for our Hck RC. 


through the inner-edge of the shell will present a lower ve¬ 
locity difference but a very high intensity due to the large 
quantity of gas integrated along the line-of-sight. Finally, a 
line-of-sight running through the outer-edge of the shell will 
display even a much lower velocity difference and a very low 
intensity due to the small quantity of gas integrated along 
the line-of-sight. In addition to these shells that cannot be 
distinguished individually on this plot, this area of the a — I 
diagram mainly corresponds to diffuse Ha emission. 


6.2 Comparison of the Ha kinematics with Hi 
results 


In this section, a comparison of our Ha kinematical results 
is done with previous Hi studies. Normally, our 5" binning 
Ha RC should give the optimal representation of the kine¬ 
matics in the inner parts of M33, without suffering from e.g. 
the beam smearing that may affect earlier Hi data. Figure [18] 
compares the Hi RC derived by|Corbelli et al.|(|2014[ l, as well 
as older Hi studies by Corbelli & Schneider ( 1997| ) an dN ew- 


ton 


(1980) and recent VLA data at 12 resolution (Gratier 


et al.||2010|| with the Ha RC obtained in this study. It is 


not surprising that the Arecibo data of |Corb clli & S chnei-| 
der (19971 suffer from beam smearing in the inner kpc and 
underestimate the rotational velocities while the 1.5 x 3.0 
resolution Newton (1980) data are quite consistent with the 
12 ,; resolution (~ 48 pc) Hi data and our 5" step (~ 20 pc) 
Ha data, at least out to 6.5 kpc. Comparing our Ha data 
to the Hi for R. > 6.5 kpc seems to suggest that we may 
overestimate the RC in the outer parts. Since the difference 
occurs at radii with the lowest number of independent bins 


(see bottom panel of Fig. 121, it could simply reflect the 
difficulty to derive accurate Ha rotation velocities by our 
fitting procedure within these regions. However it should be 


noted that the 2010 VLA data of Gratier et al. (2010) have 


lower velocities than Corbelli et al. (2014) in those same 
outer parts. 


6.3 Mass modeling results 

The extent of the spheroidal component shown in the sur¬ 
face brightness profile is very small r ^ 1.5 arcmin (^ 350 
pc); with a bulge-to-disk ratio of only ~0.04, large differ¬ 
ences between the mass models with and without a bulge 
component are not expected. In each case, a best-fit model 
and a fixed M/L model are explored. The reduced chi-square 
Xr is used to determine the goodness of the fit. 

In the pure disk case (Figure [l4| , the two DM best-fit 
models (ISO and NFW) give very similar good fits as shown 
by the values of the \r in Table [H] The redisuals (bottom 
panels) shows that the only discrepancy is in the very inner 
parts. When using the fixed M/L value of the disk, it can be 
seen that ISO gives a better fit (Xr = 1-16 vs 1.56), despite 
the fact that it still slightly underestimates the velocities for 
R < 2 kpc. This may points out that the method used to 
derive the M/L may underestimates it slightly (0.72 vs 0.83). 
This could also point out that an additional mass component 
(a small bar or nuclear disk) is missing in our models in the 
innermost parts of the disk. 

As expected, in the bulge-disk case, it is clear that the 
addition of a bulge component does not bring anything more 
than the disk only model due to its very small Re and the 
extra uncertainty coming from the M/L of the bulge. Look¬ 
ing at Figure [15] and Table [7] the conclusions are that both 
functional forms (ISO and NFW) yields very similar results. 
Because in the bulge+disk case, we add one free parame¬ 
ter (Tj, of the bulge) and because the uncertainty on that 
parameter is large since the IR color used to derive it is pol¬ 
luted by disk light with the result of underestimating the IR 
color and thus the Y;, of the bulge, we favor the results of 
the disk-only models. 

We conclude, that the pure-disk ISO models (best-fit 
and/or fixed M/L) give a slightly better representation of 
the mass distribution using the Ha rotation curve that is 
derived out to ~8 kpc with a M/L for the luminous disk 
~ 0.72 — 0.82, a core radius ~ 12 — 24 kpc and a central 
density po ~ 0.007 — 0.008 Mg .pc -3 for the dark isothermal 
halo. It is clear that when we will combine those Ha data to 
more extended Hi kinematics, we should be able to constrain 
much better the parameters of the mass distribution. 


7 SUMMARY AND CONCLUSIONS 


New Ha Fabry-Perot mapping of the nearby galaxy M33 
has been presented. The data were obtained on the 1.6 m 
telescope of the Observatoire du Mont Megantic (OMM) 
using a high resolution FP etalon (p = 765) and a very 
sensitive photon counting EMCCD camera (QE~ 90%). The 
ten fields observed cover a field of ~ 42 7 x 56 ; (10 kpcxl3.5 
kpc) with a spatial resolution ~3". The Hu regions inside 
this area are well defined spatially and spectrally. This study 
provides for the whole field, as well as for each Hu regions 
and nebulae (e.g: NGC 604, NGC 595, NGC 588, IC 137, 
IC 136) velocity and dispersion maps. The data was flux 
calibrated using Relano et al. (20131 data. 


Looking at the velocity dispersion er profile as a function 
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of radius, we found that it is essentially flat at an average 
value of 16 ± 2 km s _ . From R ~ 7.5 kpc, the velocity 
dispersions increase to ~ 20 — 25 km s -1 . This radius is the 
location of the beginning of the warp of the Hi disk where 
the twist of the position angle starts. The mean velocity 
dispersion of M33 was found to be small when compared 
to a sample of nearby star-forming galactic discs (GHASP 
sample) which have a varying from 15 to 35 km s -1 . Finally, 
a velocity dispersion versus intensity (a — I) diagram shows 
clearly a strong relation between the Ha velocity dispersion 
and the Ha intensity. 

The main aim of this study was to derive a high spa¬ 
tial resolution Ha RC for M33 in order to try to constrain 
the best functional form representing the dark matter dis¬ 
tribution between the cored ISO models and the cosmolog- 
ically motivated cuspy NFW models. The rotation center 
found is very close to the optical center. The other kinemat- 
ical parameters found are V ay s = —178 ± 3 km s -1 , PA = 
202°±4°and i = 52°±2°. The RC was derived separately for 
the approaching, receding and for both sides. The rotation 
curve was computed using a constant 5" step to perform 
the mass modeling. Besides the intrinsic error in each ring, 
we added an error to represent the asymmetry of both sides 
since this RC will be compared to axisymmetric mass mod¬ 
els. After a steep rise in the inner 1 kpc, the RC rises slowly 
to its maximum value of ~123±3 km s -1 at the last point 

(33'). 


Comparing our adopted RC to Hi rotation curves in 
the literature, we found that our Ha data agree very well 
with the Hi from the center out to ~ 6.5 kpc. The only 
exception is the first inner point of |Corbelli fc Schneider| 
(19971, which may a bit suffer from beam smearing due to 
the large Arecibo beam. With our high resolution data we 
bring more data points in the inner part of the RC, which is 
very useful for the mass modeling. On the other hand, for R 
> 6.5 kpc, this comparison suggests that our Ha data may 
overestimate the velocities in the outer parts. 

For the mass model analysis, the Spitzer/IRAC 3.6/im 
profile was used to represent the stellar disk mass con¬ 
tribution. The disk parameters found are a scale length 
Rd = 1.82±0.02 kpc and an extrapolated surface brightness 
no = 18.01 ± 0.03 mag.arcsec -2 . A bulge-disk decomposi¬ 
tion was also done even if it will not alter significantly the 
results of the mass models with a bulge to disk ratio of only 
0.04. Besides best-fit models, models with fixed M/L based 
on IR colors and population synthesis models, following the 
method used by Oh et al. (20081 and de Blok et al. (20081 
are also explored. 

For the mass models, we decided to favor the disk- 
only models because of the extra uncertainties introduced 
by adding a bulge which contributes very little (B/D = 0.04) 
to the luminous mass. In this case, the ISO models give a 
better fit despite the fact that they seem to underestimate 
the velocities in the inner 2 kpc (see the residual curves at 
the bottom of Figure 141. This may point out to the fact 
that the determination of the M/L values using the IR col¬ 
ors and population synthesis models slightly underestimates 
the M/L ratio of the disk or that an inner mass component 
has been omitted in the modeling. 

An ideal RC using our high spatial resolution Ha RC 
in the inner 6.5 kpc and a high sensitivity but low spatial 
resolution Hi RC in the outer parts should allow to con¬ 


strain much better the parameters of the mass models. This 
should be done in a subsequent study presenting new deep 
Hi observations. 
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Table Al: The Ha rotation curve of M33 
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